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Abstract 

In this paper we establish a relationship between the asymptotic form of condi- 
tional boundary crossing probabilities and first passage time densities for diffusion 
processes. We show that, under the assumption that the conditional probability 
that our diffusion {Xs,s > 0) doesn't cross an upper boundary g{-) prior to time t 
given that Xt = z behaves as (a + o{l)){g{t) — z) as z ] g(t), there exists an expres- 
sion for the first passage time density of g{-) at time t in terms of the coefficient a 
of the leading asymptotic term and the transition density of the diffusion process 
(Xs). This assumption is shown to hold true under mild conditions. We also derive 
a relationship between first passage time densities for diffusions and for their corre- 
sponding diffusion bridges. Finally, we prove that the probability of not crossing the 
boundary g{-) on the fixed time interval [0,T] is a (Gateaux) differentiable function 
of g{-) and give an explicit representation of the derivative. 
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1 Introduction 



Consider a diffusion process {Ug, s > 0) and an upper boundary g{s) (so that g{0) > Uq), 
and, for a fixed t > 0, denote by (f/f ) the process (Us) conditioned on Ut = z. The first 
main result of this paper is an explicit relationship between the asymptotic behaviour of 
the boundary crossing probabilities for this bridge diffusion process (f/|) as ^ t fl'(^) and 
the first passage density of the boundary at t by the unconstrained process. To the best 
of our knowledge, this relationship has not been observed and described in the literature. 
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We assume that our time-homogeneous diffusion process (Us) satisfies the stochastic dif- 
ferential equation 



dUs = u{Us)ds + a{Us)dWsi s > 0, (1) 

where {Ws) is the standard Brownian motion and (j{y) is continuously differentiable and 
non-zero inside the diffusion interval (that is, the smallest interval ICR such that, for 
all s > 0, t/s G / almost surely). This can be extended to some time-inhomogeneous 
processes, see Remark [7] below. We work with the transformed process Xg := F{Us), 
where 



yo 



for some Uq inside the diffusion interval of {Us)- This process satisfies the stochastic 
differential equation 

dXs = i2{Xs)ds + dWs (3) 
with fi{y) given by the composition 



see e.g. [12], p. 161. Conditions mentioned throughout refer to the transformed process 
(Xs) and its drift coefficient fi; see Remark[3]for further details on the relationship between 
results for (X^) and (Us). 

In order to establish the desired relationship between the first passage time and the asymp- 
totic conditional crossing probabilities, we rely in Theorem [T] on assumption ([U]) describing 
the asymptotic form of the boundary crossing probability for the bridge process. We then 
show in Theorem [2] that the assumption actually holds under rather mild conditions. 

These results immediately extend to first passage time densities for diffusion bridge pro- 
cesses. We show that if, for a given boundary, we know the first passage time density 
for the unconstrained process, then we immediately have the corresponding density for 
the bridge process as well (Theorem [3]). For the Brownian motion case, the results of 
Theorems [T] and [3] are closely related to those in [7] and [H] (see Remarks [T] and El) ■ 

Finally, using the methods developed when proving Theorem [21 we also show that the 
boundary crossing probability is Gateaux differentiable as a function of the boundary 
(Theorem [5]) and give an explicit representation of the derivative. 

The paper is structured as follows. In Section [2] we present the main results. Sections [3] 
toOcontain the proofs of the results. Section Hlpresenting a weak convergence result which 
may be of independent interest. Section [6] gives some examples illustrating our results in 
the Brownian motion case. 
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2 Main Results 



We denote by P^^ and E^. probabilities and expectations conditional on the process (X^) 
starting at the point Xq = x. Where no subscript is present, either conditioning is 
mentioned explicitly or the process is assumed to start at zero. We also use 1^ to denote 
the indicator of the event A. 

For the process (X^) and boundary g{s) with g(fi) > Xq, define the first passage time 

r := mi{s > : X, > g{s)} (4) 

and the transition density 



p{s,x, z) 



d_ 

dz 



Px(X, < z) 



The following theorem establishes a relationship between the asymptotic form of the 
conditional crossing probability and the density of r. 

Theorem 1. Assume that, for some < a < b < oo and G M, the boundary g{s) 
satisfies g{0) > x, 

-K'h < g{t + h) ~ g{t) < K^h, a<t<t + h<b, (5) 
and there exists a function f{t, x) continuous in t such that, for t G (a, b). 



P, 



, , snp{Xs-g{s))<0 Xt = z ] = if{t,x) + o{l)){g{t) - z) 

0<s<t 



as z 1 g{t). Then r has a density in the interval {a,b) which is given by 



(6) 



Pr{t) = 2/(^' ^)Pit' ^' 9{t)), a <t <b. 



(7) 



Remark 1. In |7] it is shown that, for a standard Brownian motion [Wt) (so that x = 
Wo = 0) and a boundary g{s) which is continuous in s G [0, t] and is left differentiable at 
t, the first passage time density at t is given by 



Prit)=bit)p{t,0,g{t)), 



where 



b{t) := lim- 



-E 



l{r>sMs)-Ws) Wt = g{t) 



sit t — S 

Thus with our Theorem [1] we have that, for the Brownian motion. 



b{t) = I hm — i P ( sup {Ws - g{s)) < 

2 z^g{t) g[t) - Z \0<s<t 



Wt = z 



Our Theorem [T] may be viewed as an alternative expression of the results of \I\ for Brow- 
nian motion and an extension thereof to general diffusion processes. 
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Remark 2. Representation ([71) can also be viewed as an extension (to the case of general 
diffusion processes and curvilinear boundaries) of the well-known Kendall's identity for 
spectrally-negative Levy processes {Xg): if g{s) = y = const, Xq = x < y and (X^) has a 
transition density, then (see e.g. pUl H] and references therein) 

Pr{t) = ^^p{t,x,y), t>0. (9) 

Remark 3. Assume that, for the original diffusion {Ug) given in ([T]), one has the asymptotic 
expression 



snp (Us -gis))<0 Ut = z ] = ifit,u) + o{l))igit) - z). 

\0<s<t / 

Then one can easily check using Theorem [T] that the density of the first passage time of 
(Us) satisfies 

Prit) = ^f{t,uy{g{t))pu{t,u,g{t)), 

where pu{s,x,z) denotes the transition density for the process (Us). 

Remark 4. We can extend the above to a class of diffusion processes with time-dependent 
drift satisfying ([6]), see Remark [7] below. 

Under mild conditions we can establish that holds for a given diffusion and boundary. 

Theorem 2. Let (Xg) be a non-explosive diffusion satisfying ([3]) with diffusion interval 
M and with jj, E such that ^ has a unique strong solution and such that there exists 
a function Q{y) satisfying 

^^\y) + ^^\y)>-Q{y), y e M, (lO) 

and 

limsup — — < — . (11) 

y^-oo y t 

Let g{s) be twice continuously differentiable for s E {a,b) , < a < b. Then there exists a 
continuous in t function f{t,x) such that, fort G {a,b), ([6]) holds as z ^ g{t). 



At this point it is worth noting that, for a linear g{s) with g{0) > x and for z < g{t), 

2 



P, ( sup (Ws - gis)) < 

KO<S<t 



Wt = z \ = 1 — exp 



t 



(g(0) - x)igit) - z) 



-(g{0) - X + oimg{t) - z) 



(12) 



as z t g{t) (see e.g. [3], pp. 64-67). This serves as a motivation for the proof of Theorem [2] 
in Section HI see Section [6] for further discussion of this example. 
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These results also extend to bridge processes, or pinned diffusions. Consider the first 
crossing time density of the processes defined as the process (X^) given in ([3]) 

conditioned to be at y at time T. Notice that, due to the Markov property, if we restrict 
our attention to the time interval [0,t], t < T, then there will be no difference between 
the distribution of the process (X^, < s < t) conditioned on being at z at time t and 
that of the process {Xs,0 < s < t) conditioned on Xt = z. That is, for any Borel set 
B C C[0,t], 



P.X Xf €5 



p. X. 



In particular, this implies that for the pinned processes the function f{t,x) is the same 
as for the unconstrained process. So, under the conditions of Theorem [2l 



P.. ( sup {Xy-g{s)) < 

yO<S<t 



XI 



P,. sup (X, - g{s)) < 

\0<s<t 

ifit,x) + oil))igit)-z). 



Xt 



Using this observation and that of Remark [7] at the end of Section O, we immediately have 
the following result. Fix y G M and denote by the density of the first hitting time 
with (Xs) replaced with (XJ) and by py{v,w,s,z) the transition density of (X^): 

p^{v, w, s, z) := P^. (Xf E dz \ XI = w) /dz, 0<v<s<T. 

Theorem 3. Under the conditions of Theorem^ the first crossing time density p'^{t) of 
the bridge process {Xf) satisfies 

Pl{t) = \f{t. a;)p^(0, X, t, g{t)), a<t<b. 

Remark 5. In [8], the authors prove an extension of the results in [Tj (see Remark [1]) to 
Brownian bridges. Denote by q^{v,w, s, z) the transition density of the Brownian bridge 
that finishes at y at time 1. Then, for (Xf ) a pinned Brownian motion and a continuously 
differentiable boundary g{s), [H] gives 

pl{t) = h{t)qy{Q,x,t,g{t)), 0<t<l, 

with h{t) defined in ([8]). Thus the result of Theorem [3] may be regarded as an alternative 
expression of the results in [8] for Brownian motion and an extension to more general 
diffusion processes. 



From Theorem [3] and the representation 

p{t,x, g{t))p{T - t, g{t),y) 



py{0,x,t,g{t)) 



p{T,x,y) 



we have the following relationship between the first crossing density of (Xg) and that of 
the corresponding bridge process. 
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Corollary 4. Assume that, for t G {a,b), (j6]) holds as z ^ 9{t)- Then the first crossing 
time densities Pr{t) and p^{t) satisfy 

y(.^ P{T -t,g{t),y) ^ , 

m = —^(j--^Pr{tl a<t<b. 

Remark 6. A result of the form holds for Bessel processes and the constant boundary 
g(t) = c > as well (see e.g. (1.1.8) in p. 429). Therefore we expect that Theo- 
rems [T], O [3] and Corollary H] will also hold for diffusions with diffusion interval (0, oo). 

Using the approach of the present paper, we can also derive an interesting result on the 
sensitivity of the boundary non-crossing probability to changes in the boundary. Denote 
by 



P{g):=Fj sup {X., - g{s)) < 

\0<s<T 

the probability that our diffusion (Xg) does not cross the boundary during the time interval 
[0,T]. It was shown in that, under broad conditions on [Xg) and g, the function P{g) 
is locally Lipschitz in the uniform norm: 

\p(g + h)-Pig)\<C{g) sup \h{t)\. 

0<t<T 

One can expect that the function P{g) will actually be (Gateaux) differentiable. The next 
theorem proves this conjecture. For simplicity, we state and prove the theorem assuming 
T = 1 (this can easily be extended to the general case). 

Theorem 5. Let (Xg) be a non-explosive diffusion satisfying with diffusion interval 
M and with fi & such that (j3]) has a unique strong solution and such that there exists 
a function Q{y) satisfying 

li\y) + li\y)>-Q{y), y g R, (13) 

and 

linisup5M<l. (14) 

y^-oo y 

Assume that g{t) and h{t), < t < 1, are twice continuously differentiable. Then there 
exists the limit 



lim e ^ 



'P{g + eh) - P{g)] = ^^^^TT^^^''^^ - r e dt) 

X Eexp {g (^-VtW^ + (7(1)) - G{g{l - t)) + ViW+g'{l) + Nt{t)} , (15) 
where the process {W^) is the Brownian meander, 

G{y) := j\{z)dz (16) 
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for some yo G M, go,u{s) := g{l — u + s),0<s + u<l, and we set 
1 



ds 



3 Proof of Theorem [T] 



Conditioning on the position of the process at time t, we have 

f9(t) 

P.(r G (t, t + h)) = F^{t e{t,t + h)\Xt = z)F,{Xt G dz) 

J ~oo 



oo 
9{t) 



P^(r e {t,t + h)\Xt = z)p{t, X, z)dz 



9{t) 



g{t)-h^/'^ J-oo 



g{t)-h^'^ rg{t) 



+0{h), 



where the last equahty follows from the proof of Theorem 3.1 in [6]. Using ([6]), the Markov 
property of the diffusion and setting A := {sup4<^<^^^(Xs — g{s)) > 0}, we have 



P.(r G {t,t + h)\Xt = z) 

= fJ sup (X, - g{s)) < Xt = z]Fj sup (X, - g{s)) > 

\0<s<t / \t<s<t+h 



Xt = z 



Jit, x){g{t) -z)+ o{g{t) - Z))P,(A \Xt = z), z< g{t). 
Define the functions 

gf{s) ■.= g{t)±K^is-t), 

so by we have 

gris) <gis) <g^{s), t<s<b, 
and hence, for < h < b — t, 

P.. ( sup (X, - g^{s)) > X, = z]< F,{A | X, = 

\t<s<t+h J 

<P,. ( sup (X, -(^-(s)) >0 
Define 

X^{z) := {g{t) - z)h-'/' ± (A + K-)h'/\ 7± W := m - z)h-'/' ± (/i - K-')h'/\ 
where 



X, = ;2 . 



fi := sup 

(;(t)-/ti/4<j/<g(f) 



fl := inf /^(y)- 

3(t)-hl/4<j^<g(t) 



P,( sup iXs-gris))>0 

t<s<t+h 



Again using derivations from the proof of Theorem 3.1 in [6] gives, as h ^ 0, 

Xt = z 

< + e2(''+^")(^(*)-^)$(x+(^)) + o{h) 

and 

Pj sup (X, - > Xt = z 

V t<s<t+h 

> Hi'iz)) + e2('^-^^)(3W-")$(7+(^)) + o{h), 

where $ denotes the standard normal distribution function, = 1 — Combining 
these results we have the bounds 

rait) 

P.(r e{t,t + h)) < / p{t, X, z) [jit, x){g{t) -z) + o{g{t) - z)) 

X (Hx'{z)) + e2(''+'^")(^(*)-")$(x+(^)) + o{h)) dz, (17) 



rgit) 

F,{r e {t,t + h)) > / p{t,x,z){f{t,x){g{t)-z) + o{g{t)-z)) 

Jg{t)-h^/i 



X ($(7- (2)) + e2(''-^^)(^(*)-^)$(7+(z)) + o(/i)) dz. 



Next we will show that 
rg(t) 



p{t,x,z){f{t,x){g{t) - z) + o{g{t) - z))<^{x'{z))dz 

g{t)-hUi 

<-hf{t,x) sup p(t,x,w) + o{h), (19) 

4 5(i)-/il/4<^<3(t) 

so that the former bound in (|T71) will yield 

Px.(r G (t,t + /i)) < ^/i/(t,x) sup p{t,x,w) + o{h). (20) 

The last claim uses the observation that the second $ term from the upper bound ( ITTj) 
admits the same upper bound as f|T9|) . since the exponential factor is 1 + o(l), and the 
o{h) term contributes a factor which itself is o{h). Using the same approach, one can 
show that the bound ffTSj) implies that 

^xir e {t,t + h))>lhf{t,x) sup p{t,x,w) + o{h). (21) 
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So we proceed to demonstrate ( fT9l) . Initially ignoring the o{g{t) — z) term, we have 

r9(t) _ 

J := p{t,x,z)f{t,x){g{t) - z)^{x~{z))dz 

rait) _ 
< f{t,x) sup p{t,x,w) {g{t) - z)^{x~{z))dz 

9{t)-h^/'^<w<g{t) Jg(t)-h}/'^ 

rait) _ 
= D / {g{t)-z)<l>{x-{z))dz, 

Jg{t)-hy^ 

where we set D := /(t, x) supg(j-)_;ji/4<^<g(j) p(t, x, w). Let y := {g{t) — z)h~^/'^ and K :- 
fi + K~ . Then 

J < ^=h / yi e-^'^'du dy 



2tT Jo \Jy-Kh^/^ 

h -h-^'^ / e-"^-"du+ - / e-"^^\u + Kh^'^fdu 



27r \ 2 Jh-^f^-Kh^/'^ 2 J_Khi/2 



D 



2V27r 

The first integral in the last line is clearly o{h) due to the term — in the exponential. 

This gives 



D 



J <^^h{l + o{l)) j w'^e-'^'^^dw + oih) 
2\/ 27r 



OO POO 



^ --hi I w^e^-^'^'dw- I w^e-^^'^'dw] +o{h). 



Again, the last integral is clearly o{h), resulting in 

J <\hf{t,x) sup p{t,x,w) + o{h). (22) 

4 a{t)-h^/^<w<g{t) 

Now consider the contribution of the o{g{t) — z) term to the left-hand side of ffT^ : 

rg{t) 



oig{t)-zMx~iz))dz. 

g(t)-h^/i 

For any £ > 0, we can take h small enough such that 



rait) _ fi 

I<e {g{t) - z)<^{x-{z))dz = e[-,h + o{h) 

Jgit)-h^/^ 



4 



using the same approach as to derive ( 122|) . so that clearly I = o{h). This proves ( fT9i) and 
hence fl2Up and fl21l) as well, as we noted earlier. Now the assertion of Theorem [T] follows 
due to the continuity of /(t, x) in t and that of x, w) in t, w. 
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Remark 7. If we consider a diffusion satisfying the stochastic differential equation 

dUt = u{t, Ut)dt + a{Ut)dWu 

then the transformed process Xt = F{Ut) with F defined in ([2]) will again have a unit 
diffusion coefficient. However, in this case the drift coefficient fi = fi{t, x) will depend on 
time. The proof of an analog of Theorem [1] in this case follows as for the homogeneous 
case, with fi and fi replaced with 

At := sup /it := inf Ks,y), 

t<s<t+h, git)-h^/^<y<g{t) t<s<t+h, g{t)-h^/'^<y<g{t) 

respectively. 



4 Proof of Theorem 



To prove Theorem [21 we first need to establish a convergence result which is similar to 
those derived in [9] and which may be of independent interest. 

In what follows, we consider processes on the time interval [0,1], that are sometimes 
pinned by their value at time 1. This can easily be changed to the interval [0,T] and 
processes pinned at time T, T > 0. 

For a > we denote by {W^) = (VT", < t < 1) the standard Brownian motion {Wq = 0) 
conditioned to arrive at a at time 1. For e > 0, set 

le{t):=e{t-l), 0<t<l, (23) 

and for < s < t < 1 define the event 



Ms,t):=\ inf (l^:-a«))>0 . (24) 

I s<u<t J 

Let Ae := v4e(0, 1) and define the conditional process iW^"^) as {W^) conditioned on the 
event A^. This will be a Markov process with transition density 

p"'"(s, y, t,z):= ¥{Wt'' G dz \ W^'^ = y)/dz 

for < s < t < 1. 

Theorem Q. As e 0, the process {Wi'^) converges weakly in the space C[0, 1] to a 
process with transition density 

p\W.z) := ^ (^1 -exp {-^}) e dz)/dz, (25) 

(l_exp{-|^}) (l-expj-f^j) / 
p\s, y, t, z) := ^ \ '-'^') 2ayx '-'^^ nW,^ G dz \ = y) dz (26) 

for z,y > and < s < t < 1. 
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Of course, the transition densities for (W^) that appear on the right-hand sides of ( l25i) 
and ( l26l) admit well-known closed form expressions (see e.g. [3], pp. 64-65). 

Remark 8. Note that the limiting process from Theorem[6]is nothing else but the Brownian 
meander on [0, 1] (see e.g. [9]) conditioned to be at a at time t = 1. This can be seen from 
comparing the transition densities for the two processes (for the transition density of the 
Brownian meander, see e.g. (1.1) in [^). 

Proof of TheoremlSi The proof follows the standard scheme based on Prokhorov's theo- 
rem. First we prove convergence of the finite dimensional distributions. Since the process 
is Markovian, it is sufficient to prove that transition densities converge. Note that if we 
take the limit as a — >■ 0, these transition densities agree with those given in Theorem 5.2 
of [H] where the process is conditioned on Wi = 0. More precisely, the following result 
holds true. 



Proof. First consider p°''^(0, 0, t, z). Using the form of the linear boundary crossing prob- 
ability for a Brownian bridge (see ( fT2l) ) as well as the Markov property of the process, we 
have for 2 > 0, as e — 0, 



Lemma 7. For < s < t < 1 and y, z > 0, 



limp"'^(0,0,t,2) =p"(0,0,t,2) 



and 



limp"''^( 



s,y,t, z) = p''{s,y,t,z). 



p'''%0,0,t,z)dz 




h(z + e(l - t)) (1 - exp {-T^(-2 + e(\ - t))a\) , 

^ ^ ^-'^ ^ '-^ (1 + o{l))F{W^^ e dz) 
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Similarly for p°''^{s, y, t, z): 
p^'%s,y,t,z)dz 

= e dz I W^'' = y)= e dz, e dy, A,)/F{W^ G dy, A, 

= F{Ae I = z, = yMWf" G dz, G dy)/F{W^ G dy, A,) 



¥{A,{s,t)\ 


\W^ = y,Wt- = z)nAit,l)\ 


= z)¥{Wt^ G dz 


= y) 




= y) 



[1 - exp {-^(y + 6(1 - s))iz + 6(1 - t))}) (1 - exp {-^(2; + £(1 - t))a}) 



1 -exp{-^a(2/ + e(l - s))} 
X ¥{Wt^ G dz I = 
(1 - exp l-^l) (1 - exp (-1^1) 
i exp I j> 



□ 



Now we will prove the tightness of the family of distributions of {W^''^). First we state 
without proof an obvious extension of Theorem 3.5 in [9j which we will use. As usual, 
denote by C[a, b] the space of continuous functions on [a, b]. Then, for s G (0, 1/2) and a 
function / G C[0, 1], denote by the restriction of the latter to [s, 1 — s]: ^/ G C[s, l — s], 

j{t) = fit)ioTte[s,l-s]. 

Lemma 8. Let {Z^, k = 1,2, . . .) be a sequence of random elements of C[0, 1]. Define the 
random elements (sZk) of C[s, l — s] as the restrictions of (Z^,) to [s,l — s], < s < 1/2. 
Then if, for any s G (0, 1/2), {sZk, k = 1,2, . . .) induces a tight family of distributions on 
C[s, 1 — s] and, for all r] > 0, 



lim lim P I sup |2'fc(t)| < rj ] =1 



and 



lim lim P sup |Zfc(t) - Zfc(l)| < = 1, 



then the sequence of distributions of [Z^, k = 1,2, . . .) in C[0, 1] is tight. 

Next we show that the conditions of Lemma [H] hold for our processes. 

Lemma 9. For any fixed a > and s G (0,1/2), the family of measures induced on 
C[s, 1 - s] by {sW^'^,e> 0) is tight. 

Proof. Denote by (VT'*') the Brownian meander on [0, 1] (see e.g. [9J for details). For sets 
B,D C C[s, l — s] and a function / G C[s, 1 — s] we use B — D and B — f to denote the 
Minkowski differences: 

B - D = {g - h : g e B ,h e D} 
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and 



B-f = {g-f:geB} 



(27) 



Define the random element (slV"^'") of C[s, 1 — s] as tlie Brownian meander (IV^) condi- 
tional on = a and restricted to [s, 1 — s]. Note that, for all y,z > and < s < 1/2, 
the joint distribution of (sW^"^'", sW^it^s) has density p"'{0,0, s,y)p"'{s,y,l — s,z). For a 
Borel set B G C[s,l — s] we also have, due to the Markov property of the process. 



P(,iy+'" e B I M^''^ = y, sWt:: = z) = P(,iy+ e B \ = v. sWt 



(28) 



Further, by comparing the transition densities for the processes and recalling definition 
fl23l). we also have 



z — es 



¥{M+ G B 



w: 



(29) 



As in [9j, for a fixed ?7 > there exists a compact set D C C[s, 1 — s] such that 

P(,iy+'" eD)>l-r]. 

The set E := (^Jg, 5 G [0, 1]) C C[s, 1 — s] is obviously also compact and so D' := D + E 
is compact, too. Clearly, for any e G (0, 1), := D + sh C D'. Then we have, as £ J. 0, 
setting := {—e{l — s), oo) x {—es, oo), / := (0, oo) x (0, oo) and using (l29l) . 

F{sW'''' G D') > F{sW'''' G D,) 

= I p(,iv"'" G De I ly;-" = y, PFfi", = z)^{w^'' G H^f:", g rfz) 

= y ¥{sW^ eD,-,l,\Wt = y + e{l-s),W+^^ = z + es) 

X p(iy;'" G dy, vrfi", G dz) 

= y p(,iy+ G D I w^/ = y', = z') 

X P(iy;'" G dy' - e(l - s), Wfl", G rf^' - es) 

jF{sW+ eD\w+ = y', iy+„, = /)p(iy+'" G dy', G d^') 
= j F{sW^''' G D I = w^it:" = z')F{w^''' G rfy', w^it:^ g rf/) 

=P(,Ty+'" G > 1 - r], 

where the convergence is justified by Scheffe's theorem (see e.g. Theorem 16.12 in [2], 
p. 215) and the second last equality uses Thus there exists an such that F{sW"''^ G 
D') > 1 — 27] for e < En which completes the proof of the lemma. □ 
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Lemma 10. For any 7] > 0, 



limlimP ( sup \Wt''\ < t] ] = 1 (30) 

0<i<s 



and 

limlimpf sup \Wt''-a\<ri]=l. (31) 

s^Oe^O \l-s<t<l J 

Proof. The proof uses an argument similar to the one demonstrating Lemma 5.4 in [U]. 
To estabhsh (IHUj) . first note that 

P ( inf Wt'' > -e] =1, 

and so we just need to consider 

P ( sup Wt'' <v] = P ( sup W^^ <ri Ae] = f( sup W^^ < rj; aA /p(v4,). 

\0<t<s J \0<t<s J \0<t<s J I 

Denote by q{t,y,z) = ^¥y(Wt < z) the transition density for the Brownian motion 
process. By conditioning on the value of the process at time s, using the Markov property 
and the known closed form expressions for both the joint distribution of the maximum 
and minimum of the Brownian bridge (see e.g. (1.15.8) of |3j) and the distribution of the 
maximum of the Brownian bridge given in ( fT2l) . we obtain 

P( sup < 7]] A, 

\ 0<t<s 

= r P ( sup W^^ < 7]] A,{0, s) = z) ¥{A,{s, 1) | = z)F{W^ G dz) 

J~e{l-s) \0<t<s ) 

r P ( sup W^^ < T], inf > -e{l - s) = z) P(y4,(s, 1) I = z) 

J-e{l-s) \0<t<s 0<t<s J 

X p(iy; G dz) 

V \q(s, 0,z + 2k[r] + e{l - s)]) -q(s,0,z + 2k{r] + e{l - s)) + 2e{l - s] 

-e(i-s) L V y V 

(^1 - exp {-^^(^ + ^(1 - ^))}) {q{s, 0, z)y'F{W: e dz). 



> 



'-e{l-s) 
X 

Now dividing both sides by 



¥{Ae) = 1 - exp{-2ea} 

(cf. ( IT^ ) and using L'Hopital's rule, the dominated convergence theorem and the fact 
that ~ N{as, s(l - s)), we obtain 



limpf sup Wt' <v] > r V - — ^(2r/fc + z)e-2^^('''=+")/^ 
^-0 \o<t<s J Jo as 



X ( 1 — exp \ \ ) g(s(l — s), as, z)dz. 
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As in the proof of Lemma 5.4 in [9], the integral of the term with k = approaches 1 as 
s — i> 0. The contribution of the terms for k > 1 is clearly positive, so to complete the 
proof it is sufficient to show that the absolute value of the contribution from the terms 
with k < tends to as s ^ 0. Setting m = —k, the absolute value for these terms is 
bounded above by the expression 

2r]{l -s) n f r 2za 



as 



(^1 -expj-^^l^ g(s(l-s),as,z)^me-2^'"('""-")/^d2. (32) 



m=l 

The contribution of the term with m = 1 is bounded above by 
27]{l-s) rv 



as 



q{s{l- s),as,z)e-^^''^''-'^/'dz 







= 2^(l-^) ,-2.(.-a) U _v0^2s)+as_\ _ ^ / 2r^{l-s)+as j j ^ 
as \ \ y^^il^ J \ y^Jil^ 

as s ^ 0. The same approach shows that the m = 2 term also gives a contribution which 
is o(l) as s — > 0. For the remaining terms, observe that 

oo oo 
m=3 m=3 

We may assume without loss of generality that s is small enough such that the summand 
in the above expression is a strictly decreasing function of m for m > 3. We then have 
that 



m=3 



where the last relation uses Mill's inequality. Thus the absolute value of the sum of all 
terms for m > 3 in ( !32l) is bounded above by 

~ r (l- exp qisil - as, z)dz < ii^e'^^^/^ = o(l) 



•qa Jo \ I 1 ~ J / Va 

as s — > 0. This completes the proof of ( l30l) . 

Now consider fl3T]). We have 



P( sup \Wt'' -a\ < 7] 



l~s<t<l 



P ( inf Wr -a>-ri, sup - a < r] 

l-s<t<l l^s<t<l 

P( inf Wt^>-ri + a, sup 1^° < 77 + a; A, ) /P(A,) =: Ji/ J2. (33) 

\l-s<t<l l-<i<t<l 
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For Jl, conditioning on the value of Wi_g gives (without loss of generality, we take a—r] > 
— if this were not the case, we would replace the lower limit of integration with 0) 



Jl 



a+rj 



a—rj 



P inf Wt> -r] + a, sup Wt < i] + a; A, 

l-s<t<l l-s<t<l 



a+ri 



a—Tj 



P ( inf Wt> -Tj + a, sup Wt < r] + a 

l-s<t<l i_s<t<l 



Wis = y,Wi = a 

X ¥{Wi-s edy\Wi = a) 
Wi^s = y,Wi = a 
X P (A,(0, l-s)\ Wis = y) nWi-s edy\Wi = a) 



Since P (A,(0, 1 - s) | Wi_, = y) = 1 - exp {-26{y + e(l - s))/(l - s)} and J2 = 1 
exp {—2ea} (see (fT2ll ). using the dominated convergence theorem in (!33l) gives 

limpf sup iVFt'''" - a| < r/^ = /" V(H^i_, e I W^i = a) , ^ , 

e^O \l-s<t<l J Ja~v a(l-s) 



X P ( inf Wt>-ri + a, sup IVt < + a | Wis = y,Wi = a] . 

l~s<t<l i^s<t<l 



This clearly tends to 1 as s — > 0, completing the proof of Lemma [TOl 



□ 



Combining Lemmas [814T0] completes the proof of Theorem El 



□ 



We now prove Theorem O 



Proof of Theorem 0. The proof of Theorem [2] can be split into three steps. In the first 
step, we change the measure to express the conditional probability of the event on the 
left-hand side of ([6]), 

A:=(sup(X,-(^(s))<ol, (34) 
Lo<s<t J 

in terms of the expectation of a (nice) functional of the Brownian bridge process (the 
Brownian motion starting at x at time and pinned at z at time t), using the approach 
employed in [1| and later in [6j. Now the underlying "pre-conditional" process (the Brow- 
nian motion) is space-homogeneous, and therefore we can "straighten" the boundary g{s) 
by switching from the canonical process (X^) to the process {Xg — g{s)) (the second step). 
In the transformed space, the original set A becomes 

i := I sup X, < I , (35) 

lo<s<t J 

and the conditional process can be thought of as obtained from the diffusion (X^) with 

dXg = -g\s)ds + dWs (36) 
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{(Ws) being a Brownian motion), Xq = x := x — g{0) < 0, conditioned on Xt = z := 
z — g{t) < 0. We again change measure (now in the transformed space) to express the 
desired expectation in terms of another one, again for a Brownian bridge process — now 
starting at x at time and finishing at z at time t — over the event A. In the third step, we 
re-write the latter expectation as the product of the (known) probabihty of the Brownian 
bridge to stay below zero (this factor will have the desired behaviour az{l+o{l)) as z ^ 0) 
and the conditional expectation, where the conditioning now includes the event A as well. 
It remains to observe that the Brownian bridge is 'time-reversible' (if {Yg, < s < t) is 
a Brownian bridge 'pinned' at times s = and s = t, then {Yt_s, < s < t) is also a 
Brownian bridge pinned at these same times), and so the behaviour of the last conditional 
expectation as 5 j can be found using our Theorem [6l 

Now we will make the above-outlined argument more precise. 

Step 1: Let Q^, denote the law of the Brownian motion (Wg) in C[0,t] with Wq = x and 
the law of (Ws) with Wq = x and Wt = z. We begin by recalling the following result 
from [T]. 

Lemma 11. LetF^ denote the law of{Xs) governed by ([3]) starting at Xq = x and pinned 
by Xf = z. Then, for any B E cr [X^ '■ u <t), 

p[t,x,z) '- -' 

where q{t,x,z) is the transition density for the Brownian motion and E^. denotes expec- 
tation with respect to the probability Q^, G{y) is defined in f|T6l) and 

1 /■* 

N{t) :=-- if-i'iX^) + /i'(X„)) du. (37) 
^ Jo 

Note that the notation of the above lemma tacitly assumes that (Xs) is a canonical process 
on the sample space C[0,t], and we continue with this assumption throughout the proof 
of Theorem [21 

Hence for the event A defined in flMl) we have 



P. ( sup (X, - g{s)) <0 Xt = z)= F^iA) = ^i^lMe«(^)-G(.)|,. ^^mt)^^-^ (gg^ 

\0<s<t J P{t,X,Z) 

Note that the changes of measure used here and in the sequel using Girsanov's theorem are 
justified since under these measures (X^) is non-explosive and fi{y) (and g{s)) is locally 
bounded, implying the conditions of Theorem 7.19 of [H], as well as Assumptions (I) an 
(II) of Theorem 7.18 (with our assumption of the existence of a unique strong solution to 
(ED), hold. 

Step 2: We now transform the space, defining ip : C[0,t] C'[0,t] by ipif) = f — g- This 
mapping induces a new measure on the space: for a Borel set B C C[0,t], 

Q,m = QAB + g), 
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again using the Minkowski difference notation defined in ( 127|) and x = x — g{0)- In the 
same way ip also induces the measure Q| (with expectation E|) from Q^. This gives 



where 



El [e^Wl^] =1:! [e^Wl^' 



[fi'{Xs + g{s))+fi\Xs + g{s))] ds 



and A is given in ( l35l) . Note that under the measure Q| the canonical process (X^) is no 
longer a Brownian motion. However, as stated previously, we can think of the conditional 
process as being obtained from the diffusion (X^) defined in (1361) . We change measure 
again such that under this new measure, Q^, the process (X^) is again a Brownian motion 
starting at Xq = x. By Girsanov's theorem. 



Cf.-- 



exp <^ - / g'{s)dXs 



{g'{s)rds 



Using Ito's formula we have 

f g\s)dX, = Xtg'{t)- fg"{s)XJs. 
Jo Jo 



This gives 



Ct = exp <j -Xtg'it) + I g'\s)Xsds - ^ / {g'isjfds 



Let (with corresponding expectation E^.) denote the law conditioned on Xt = z; 
this is clearly the distribution of the Brownian bridge on [0, t] pinned at x at time s = 
and at z at time s = t. Applying the same reasoning that leads to the assertion of 
Lemma [TT] (see [Ij) then gives 



q{t,x,z) 



E 



p{0,X^t,z) ^_-^,^^)-z 

q{t,x,z) 



(39) 



where p(0, x, t, z) is the transition density of the (time-inhomogeneous) process (X^) under 
Q^, and we set 



Nit) ■■=-\l f^' {X, + g{s)) + /i^ (X, + g{s)) 



ds 



Combining fl38|) and fl39l) then gives 



P,. sup (X, -g{s))<0 

\0<s<t 



Xf = z 



+ j\\s)Xds-\j\g'{s)fds. 



q(t, X, z) p{0, X, t, z) ^G(z)-G(x)-zg'(t)'^^ 

p{t,x,z) q{t,x,z) 
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Step 3: Now condition on the event A. From ffT^ . we have Qs(^) 
which gives 



1 — exp{~2xz/t}, 



P,( sup (X, - < 

0<s<t 



Xt 



q{t,X,z) p{0,X,t,z) ^G(z)-G(x) 

p{t,x,z) q{t,x,z) 

q{t, X, Z) P{0, X, t, ~Z) ^G(z]-G(x)~z 

p{t,x,z) q{t,x,z) 



e"'~' ""^'^^^ ( 1 - exp <^ 



t 



A 



(1 + 0(1)) 



Nit) 



A 



(40) 



as z ^ g{t) (i.e. z ^ 0). Due to the above mentioned time-reversal and symmetry properties 
of the Brownian bridge, we obtain that 



.Nit) 



A 



Mi) 



A 



(41) 



where 



Kit) = -ll + 9{t - s))+AXs + git - s))] ds 

+ f g"{t - s)XJs - \ !\g\t - s)fds. (42) 
^0 Jo 

Thus we are now in the situation of Theorem [6] (with an appropriate change of scale), 
conditioning on our event A being equivalent to conditioning on the event ^^(0, 1) (see 
f l24p ) in the theorem (note that the initial value is now z ^ 0). To complete the proof, we 
need to show that the expectation (14 ip converges to a finite limit as i | 0. For a fixed 
large H > 0, let D := {info<s<tXs < —H} and, as usual, denote the complement 
event. We have 



e: 





A 


= n 




A 






A 



(43) 



In the integrand in the second term on the right-hand side of (l43l) . iV(t) is clearly a 
bounded continuous function of X. (in the uniform topology) on A. Since D has zero 
boundary under the limiting distribution, we see, using the weak convergence result of 
Theorem [Ul that this term converges to a finite limit as z t 0. We complete the proof 
of Theorem [2] by showing that, as H —>■ oo, the first term on the right-hand side of fHHl) 
converges to (which is basically equivalent to uniform integrability of e— under the 

distributions QI{-\A), z e (-1,0)). 



We begin by considering the distribution of the minimum of (X^) for s G [0,t], under the 
measure and conditional on A. Define 7 := inf{s > : Xg = x}. Under Q^, we clearly 
have that 7 < t a.s. Then for y < x, observe that 



inf Xs <y 

0<s<t 



A 



inf Xs <y 

0<s<t 



A; 7 



u 



(44) 
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Using the strong Markov property of (Xg) and ffT^ (recalling that, under Q^, the process 
[Xg) is a pinned Brownian motion), we have that 



inf Xs <y 

0<s<t 



A;'y = u 



inf Xs <y 

u<s<t 



A;'y = u 



inf Xs < y, sup < 

u<s<t u<s<t 



x„ 



X 



< 



inf Xs<y 

u<s<t 



Xu = X 



sup Xs < 

u<s<t 



sup Xs < 

u<s<t 



Xu = X 



x„ 



X 



_ exp{-^{y-xy} ^ e-2fa-i)Vt 
1 — exp{ — ^x^} — 1 — e-2*V* 
Hence an upper bound for fl44l) is given by 



inf Xs <y 

0<s<t 



A < 



:(7 e du\A) 



_ p-2£2/t ■ 



(45) 



1 _ g-2£2/t 

Now consider iV(t) in the first term on the right-hand side of (l43l) . Set 'g{t) := maxo<s<t (?(s) 
and g{t) := mino<s<t 5^(5). Using (H2ll . assumption (fTOi) and the assumption that (7(5) is 
twice continuously different iable, we see that, on the event A, 

N{t)<-c mi Xs-l inf [^'(y) + ^^(y)] 

2 mio<s<tX,+g{t)<y<g{t) ^ 



0<s<t 



t 



< -c inf Xs + -Q inf Xs + g{t) 



0<s<t 



Q<s<t 



for some c > 0, where we have assumed, without loss of generality, that Q{y) is a de- 
creasing function for y < 0. Thus using ( |TT|) . for large enough H and some r < 2/t, we 
have 

X(t) < r ( inf XJ on D. 

\0<s<t 



Hence, using ( l45l) . 



A 



< 



-H 



inf Xs <y 

0<s<t 



inf Xs<y 

0<s<t 



1 -H 



-H 



2rye 



ry 



inf Xs < V ] dy 

0<s<t " - ^ ' ^ 



< 



ye 



rj/2-2(j/-i)2/t^. 



(46) 



which clearly vanishes as H ^ —00. Thus the expectation in ( HOi) converges to a finite 
limit and hence we have 

( sup (Xs - g{s)) < Xt = z]= (fit, x) + o(l))5 = (/(t, x) + o{l))ig{t) - z), 

\0<s<t / 

completing the proof of Theorem [21 



□ 
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5 Proof of Theorem 
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The proof of this theorem can be divided into four steps. For the first three steps, we 
assume h{t) > 0, < t < 1. In the first step, we observe that the difference P{g + eh) — 
P{g) can be written as an integral by conditioning on the first crossing time r of g. In 
the second step we follow a similar scheme to the proof of Theorem [21 transforming the 
integrand so it is written as the product of an expectation of a functional of the Brownian 
meander and a well-known boundary non-crossing probability for the Brownian motion. 
In the third step, we calculate the limit of the ratio of the thus obtained expression to 
e as e — >■ 0. This involves careful treatment near the right end point of the integration 
interval. Finally, in the fourth step, we show how to extend the results to general h which 
are twice continuously differentiable. 

Step 1: Clearly, the difference P{g + eh) — P{g) is the probability of (Xg) crossing g at 
some time prior to time 1 without ever crossing g + eh on [0, 1] (as we said, we assume 
here that h > 0). We condition on the time remaining until time 1 after the first crossing 
time of g and use the strong Markov property to obtain that, for e G (0, 1), 

P{g + eh) - P{g) = [ F,{1 - r e dt)Fj sup (X, - g{s) - eh{s)) < Xi^t=g{l-t)^ 
Jo Vi-*<^*<i / 

= [ P,(l -re rft)P^(i_t) ( sup (X, -g{l-t + v)- eh{l -t + v))<o] 

Jo \0<v<t / 

= / + / =: Ji + J2. (47) 



' e 



Step 2: We will show that Ji = o{e) in Step 3. For J2, we continue in a similar manner 
to the proof of Theorem [21 Note, however, that in this case we do not condition on the 
end point of the process. Define 

ge,t{v) := g{l-t + v)+ eh{l -t + v), 0<v<t, 

and set 



A{t) := <^ sup (X, - ge,tiv)) < 0) 

\^0<v<t 

Using Girsanov's theorem to change to the Brownian motion measure, we have 

where ]Eg(i_t) denotes expectation with respect to the measure Qg(i-t) defined in the proof 
of Theorem [2] and N{t) is given by ( |371) (see the proof of Theorem [2] for the justification 
of the use of Girsanov's theorem and Ito's lemma). Again we transform the space using 
the function ip^^t : C[0,t] C[0,t] defined by ipeAf) ~ f ~9e,t- This induces the measure 
Qs(t) with expectation K^i^t), where 6{t) := —eh{l — t), such that 

F ^ \pGiXt)~Gig{l-t))+N(t)^ 1 _ ^ \ GiXt+g,At))-Gigil-t))+Ne,t(0,t)-a _ 
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where 



and 



[^'{Xy + ge,u{v)) + ii^{Xy + ge,u{v))\ dv 



A{t) := { sup X, < ^ . 

Again we can think of {Xy) under this new measure as driven by the stochastic differential 
equation 

dXy = —g'^ ti'^)dv + dWv, 



for a Brownian motion (Wy). Changing to the measure Qs(t) (with corresponding expec- 
tation ^s(t)) such that {Xy) is again the Brownian motion gives 



P,(i_i) (Ait)) = Es^t) 



G{Xt+9e.t{t))-G{gil-t))-Xt9',^tit)+Ne.t{0,t) 11 - 



where 



iV,,„(s,t) :=iV,,„(s,t)+ / g':jv)Xydv-^ j {g'.^Mfdv. 



We then condition on which yields 



^G(Xt+ge,tm-G{9(l-t))~Xtg',Ai)+N,^t{0,t) 



A{t) Qsit)[A{t)). (48) 



The last factor is now the probability for the Brownian motion to stay below a fixed level 
on [0,t]. Hence, uniformly in t G (£'^^^, 1), we have 



Q,(,)(i(t)) = 2$ (^^^^1-1 
as e — >• 0, see e.g. (1.1.4) in [3], p. 153. 




(49) 



Step 3: We now divide the right-hand side of fH7|) by e and take the limit as e — > 0. 
Let c denote an upper bound for the density '^xi)- — t ^ dt)/dt on [0, 1] (which can be 
obtained, for example, using Theorem 3.1 of [6j). Then we clearly have 



Ji < c / dt = ce^l'^ = o{e) 
Jo 



Using (jlH]) and 1^, we see that 



J2 

e 



TT J 1,3/2 



h{l-t)E, 



,G{Xt+geAt))-G{g{l-t))-Xtg',Jt)+N,,t{0,t) 



Ait) 



X ^ ^(1 + 0(1))- (50) 
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Using the weak convergence result of Theorem 2.1 in [9], together with assumption ( fT3i) 
and a similar uniform integrability argument to the one used in the proof of Theorem [2] 
(see (H3|) - (l46|) : that the right-hand side in condition (fT^ differs from that in (ITT!) is due 
to our dealing with the maximum of the Brownian motion, while in Theorem [2] we dealt 
with that of the Brownian bridge — which has a thinner distribution tail), we have that, 
as e — > 0, the expectation in ( !50l) converges to the expectation under which the process 
{—Xs) is the Brownian meander on [0,t] (that is, the Brownian motion (1^5,0 < s < t) 
conditioned to remain positive on [0,t]). If we denote this process by {W^'^'^), then, by 
the Brownian scaling, we have 

Thus the expectation on the right-hand side of ( l50l) converges to the expectation in the 
statement of Theorem El which completes the proof of Theorem for h >0. 

Step 4: Now consider a general h{t) satisfying the conditions of the theorem. Using the 
standard notation y~ = — min{?/, 0} for the negative part of ?/, we have 

P{g + sK) - P{g) = {P{g + eh) - P{g - eh-)) - {P{g) - P{g - eh')) . 

Each of the two terms on the right-hand side can then be evaluated as per Steps 1-3, 
the "lower" of the two boundaries (originally it was g due to the assumption h > 0) now 
being g — eh~ . Next we observe that the main terms in the respective expressions for 
J2/e are continuous functions of e, the limits of them having the form of the right-hand 
side of (1151) with h{t) replaced with max{/z.(t), 0} and max{— /i(t), 0}. This completes the 
proof of Theorem 



6 Examples 



Example 1: Brownian Motion and a Linear Boundary. In this example we will illustrate 
the results of Theorems [T] and in the case of the Brownian motion (Wg) when g{s) 
is linear, with g{0) > x. Recall that (fT2l) shows that the conditions of Theorem [T] are 
satisfied with f{t, x) = j{g{0) — x). Thus the theorem implies that 



Pr{t) = ^j{g{0)-x)p{t,x,g{t)) = 1((^(0) -x)- 



,-(9(t)-x)V(2t) 



which is a well known result, being a special case of Kendall's formula 
of [3], p. 250. 



see e.g. (1.1.4) 



Next assume that g{t) = ai + bit and h{t) = 02 + 62^, < t < 1, where bi, 02, &2 £ and, 
without loss of generality, ai > 0. Then the assertion of Theorem [5] reduces to 



lim e 



P{g + eh)-P{g) 



2 as + 62(1 -t) 



P^(l -re dt)E 



(51) 
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On the other hand, we know (see e.g. (1.1.4) on p. 250 in [3]) that, for the boundary 
a + bt, a,b > 0, 

P(-cx), a + bt) = $ (-a - b) + e-^^^^ (-a + b) , 

(52) 

Pr{t) =^e.p{-i^}, 
while, from e.g. (1.1) in [9], we have that 

F{W^ e dy) = ye~y'/^dy, y > 0. (53) 



In the special case 6i = 62 = 0, we can evaluate both sides of ( 15T1) and confirm they each 
give 

-ai/2 



— 026 
TT 



In the general case, we use f loHj) to evaluate the required Laplace transform: 

E[e^<] = 1 + v^Ae^'/'$(-A), XeR. (54) 
Using (!52|) and fl54|) . in this case (15T!) is equivalent to the relation 



V vr 



ai /"^ a2 + 62(1 -t) / (ai + &i(l -t))^ 



02 



TT 7o V^(l - t)3/2 ^ i 2(1 - t) 2 

X (^1 + V2^6ie*'''/^$(-v^6i)) dt. (55) 

It is not immediately clear that this identity holds true. However, for given values of ai, 
02, bi and 621 we can numerically integrate the right-hand side of (1551) to confirm (1511) 
holds for these values. For example, using the values Oi = 02 = 61 = 62 = 1, both sides of 
(l55l) give 0.379. Alternatively, for the values ai = 1, 02 = —0.5, bi = —1, 62 = 2 (in which 
case the "increment" h{t) assumes values of both signs on [0, 1]), both sides of (1551) give 
0.442. 

Example 2: Brownian Motion and a Daniels Boundary. A general Daniels boundary is 
given by 



9d{s) := ^ - ^log (^y + + n,e-^''/^^ , (56) 

where 5 7^ 0, ki > and K2 G M are subject to n\ + 4/t2 > 0. For this boundary, the first 
crossing time density is given by 
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(see [5] for further information). Theorem [T] therefore gives 



(57) 



Figure [T] shows estimated values of the probabihty 



( sup {Ws - gois)) < 

^0<s<t 



(5^ 



using simulation, as a function of goit) — z. We also estimate the value of f(t, x) predicted 
by this simulation (using only the simulated values at the points < goit) — z < 0.1) 
and plot the corresponding linear function for comparison. The parameter values used 
are 6 = ki = K2 = 0.5, t = 1 and a; = 0. Simulation of the Brownian bridge process 
was performed using the fact that if (Bg) is a Brownian motion conditioned on Bf = z, 
Bq = X, then 

£ X + - {{z - x) - Wt) + Ws, 0<s< t, 

see e.g. [Hj, pp. 64. We used 10^ simulations for each value of goit) — z; the step size 
is 10~^ for < t < 0.99 and 10~^ for 0.99 < t < 1 (we use a smaller step size closer 
to t = 1 as this region is prone to the most simulation error). The regression gives the 




0.05 0.1 0.15 0.2 



gD(t) -z 



Figure 1: Simulated values of fl3^ as a function of g£,{t) — z for the Daniels boundary 
fl3U|) with 6 = Ki = K,2 = 0.5 using x = and t = 1 (solid line). The dotted line shows the 
linear function corresponding to the regression of the values for < goit) — z < 0.1. 

estimated value of /(1,0) as 1.35, in comparison ^7\\ gives /(1,0) ~ 1.33. The slightly 
higher value from simulations is to be expected since simulating continuous stochastic 
processes using discrete simulation points over-estimates the probability the process will 
stay below a given boundary, as the simulated process is unable to cross the boundary 
between two consecutive time steps. 

Note that the results of Corollary H] immediately give the form of the first passage time 
density for the conditional process in this example. 
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